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Abstract 

We solve two-dimensional large- QCD in the presence of a nonzero baryon number B, and for 
arbitrary quark mass m and volume L. We fully treat the dynamics of the gluonic zero modes and 
check how this affects results from previous studies of the B = and B = 1 systems. For a hnite 
density of baryons, and for any m > 0, we find that the ground state contains a baryon crystal 
with expectation values for {tpj^tp) that have a helix-like spatial structure. We study how these 
evolve with B and see that the volume integral of {tpil^) strongly changes with the baryon density. 
We compare this emerging crystal structure with the sine-Gordon crystal, which is expected to 
be a good approximation for light quarks, and hnd that it is a very good approximation for 
surprisingly heavy quarks. We also calculate the way the ground state energy E changes as a 
function of the baryon number B, and find that for sufficiently large densities the function E{B) 
is well described by the equation of state for free massless quarks, thus suggesting a quark-Hadron 
continuity. From dE{B) / dB we calculate the quark chemical potential as a function of B and 
see that the baryons repel each other. The way /i depends on B also allows us to translate our 
findings to the grand-canonical ensemble. The resulting phase structure along the ^u-axis contains 
a phase transition that occurs at a value of ^ equal to the baryon mass divided in A^, and that 
separates a /i-independent phase with intact translation symmetry from a /i-dependent phase with 
spontaneously broken translation symmetry. Finally, our calculations confirm the presence of a 
partial large- Eguchi-Kawai volume independence, as described in Phys. Rev. D79: 105021, that 
arises only if one treats the gluonic zero modes correctly. 
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I. INTRODUCTION 



A first principles calculation of QCD properties in dense environments is a challenging 
problem facing the community working on non-perturbative aspects of QCD. This is mostly 
due to the severe sign problem inflicting lattice simulations at these densities. Trying to avoid 
the sign problem typically leads one to study relatively low densities Other approaches 
such as analytic weak/ strong coupling expansions come with their own shortcomings : weak 
coupling expansions rely on the smallness of the effective interaction between he quasi- 
particles of dense QCD, but these interactions really become weak only for asymptotically 
large chemical potentials j^. On the other hand, lattice strong-coupling methods 4] are, of 
course, far from the continuum limit. Furthermore, the approach of using the gauge/gravity 
duality is also not free of difficulties [5|] : it is truly systematically controlled only when the 
curvature scale associated with the classical gravity solution is small in its natural units. 
Unfortunately, the classical solution for a single or few baryons has strong curvatures and 
can get 0(1) corrections from higher derivative terms in the effective low energy action, 
which are currently unknown. For a large number of baryons one also needs to assume small 
curvatures, and this restricts one to a translational invariant states which may have little to 
do with the true ground state. 

Given the above theoretical status, we find it surprising that the soluble system of large- 
QCD in 1 + 1 dimensions - the 't Hooft model gI - has not been generally solved for such 
densities. For example, the interesting works of Refs. [3] either restrict to zero quark mass, 
or expand in the value of the m. Unfortunately, for such values of m the two-dimensional 
baryon is almost massless and may behave very differently from a four-dimensional baryon. 
A different theoretical approach to this problem is presented in Ref. 8| and invokes the 

u n 

'Eguchi-Kawai' equivalence [9(. We defer a short discussion on Ref. [8| to our summary. 

The aim of this paper is therefore to solve two-dimensional large-A^ QCD for general 
values of the quark mass m, the baryon number B, and the spatial volume L. 

The outline of this paper is as follows. Section [III describes the results of our previous 
publication [10] where we derived an IR and UV regularized form of the 't Hooft classical 
Hamiltonian, and its corresponding 't Hooft equations. Besides depending on m, on B, 
and on L, this construction depends on the lattice spacing a and on an additional technical 
parameter which we denote by M and that needs to be taken to oo. The role of M is 
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to discretize a certain continuous variable on which the solution of the 't Hooft equations 
depends. Comparing the construction of Ref. jlO| to other works such as Ref. [ill] we see 
that these instead have set M = 1. Such choice is inconsistent with the dynamics of the 



10|). Thus, in 



spatial glue fields, especially when the volume is small (see [12|], [13[, and 
Section UTI] we check what is the effect of increasing M from 1 on the properties of the ground 
state for 5 = and B = 1, which were already calculated previously in Ref. [ll^. Section HVl 
contains our main results and describes the way a crystal is formed if the system contains a 
variety of baryon numbers in the range < -B < 30. In Section |V] we confirm the presence 
of a partial volume reduction anticipated in Ref. [10|]. Section |VT] summarizes our results. 



II. THE CLASSICAL 'T HOOFT HAMILTONIAN Hp AND ITS MINIMIZATION 

In Ref. [l^ we made a first step towards solving large- QCD in 1 + 1 dimensions 
for arbitrary values of baryon number, volume, and quark mass. For the purpose of self- 
completeness we use the current section to repeat the main results of that paper. 

We derived a regularized form of the 't Hooft classical Hamiltonian which is both UV 
and IR finite by using the canonical formalism defined on a one-dimensional spatial lattice 
of spacing a and volume L = aLs (here Lg denotes the number of lattice sites). The fermion 
discretization of our choice was a single species of staggered fermions, which gives rise to a 
single Dirac fermion in the continuum (i.e. there is no residual doubling).^ 

Since our system is IR regularized by a finite box, we cannot gauge away the spatial glu- 
ons completely. Specifically, a set of — 1 zero modes which correspond to the eigenvalues 
of the spatial 'Polyakov loop' remain dynamical. Taking this into account and using the 
coherent state formalism of Ref. [l^ we obtained a form for the 't Hooft Hamiltonian that 
depends on the expectation values of the Polyakov loops and of the fermion bilinear opera- 
tors. Minimizing this classical Hamiltonian in a 1/A^ expansion we saw that the Polyakov 
loops have zero expectation values and that the classical Hamiltonian becomes a functional 
of an infinite set of basis functions. Each such basis contains Lg orthogonal scalar functions 
that span the L^— dimensional space furnished by the spatial coordinate x G [1,-^^^]. We 

^ This regularization is self-consistent only if Lg is chosen to be even. 
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denote these functions by (pa{x): 

a = l,2,...,oo ; n,xe[l,L,], (2.1) 

with 

E(C(^))^C(^) = Sr^m, (2.2) 

X 

Ec(^) my)r = 5xy. (2.3) 

X 

In terms of the functions the regularized 't Hooft Hamiltonian becomes (here we also 
regularize the infinite number of functions to be the finite number M) 



HFi(t>)lN = lim 



Ti^T {(4pL« + c.c.)+™(-im4 



E E 



4^«M2t^^.^^^^ 4sin^ g(^ + 0/2 



\{B + Ls/2) 1 1 
/I r /\/r2 2^ 



i^.M2f^4sin2 (g(^^ + /)/2 



(2.4) 



Here and below the bare mass m and bare 't Hooft coupling A are dimensionless and related 
to the corresponding dimensional quantities m, A via 

A = a^A, (2.5) 
m = am. (2.6) 

Also, the set of M 'density matrices' p°;.y is given by 

B+Ls/2 

Ply^ E €{x)€\y)- (2.7) 

n=l 

The 5-dependence of Eq. (12.71) . together with the way Hf{4>) depends on B, is the way 
the baryon number enters the discussion (in the Hamiltonian approach one constraints the 
quantum Hilbert space to obey Gauss law and for the SU (N) gauge theory that we study 
this includes fixing the global U{1) baryon charge B).^ Finally, note that the prime on the 
sums means that the terms with a = b and / = Lg are excluded (/ generally gets integer 



^ In our numerical studies we find that, while irrelevant for the minimization oi Tip, the last term in 
Eq. (|2.4p is crucial to include in order to get the correct baryon mass. 



3 



values between 1 and ) . This exclusion is a direct result of our IR regularization and of the 
correct treatment of the zero modes. It is exactly this exclusion that becomes the principle 
value prescription often used in differently regularized treatments of the 't Hooft model.^ 

Let us make an important remark already alluded to in the previous section: in [lO(] 
we showed that the choice of M = 1 (instead of M = oo), that was made elsewhere (for 
example in [ui), is effectively equivalent to incorrectly treating the spatial gluonic zero 
modes and that it amounts to choosing the expectation values of the spatial Polyakov loops 
to be nonzero. This is wrong and leads to erroneous results. In particular, it implies that 
the volume dependence of the Hamiltonian is in contradiction with general arguments of 

n I 

the sort of [9[ (for a detailed discussion on this point see p^]). In the current paper we test 
how increasing M from 1 to larger values of 0(15 — 25) fixes this problem, and we do so for 



different values of B. This extends the 5 = discussion of 



Ref. [u] for 5 = 1. 

Going back to the derivation of the 't Hooft equations, we perform the variation of Tip 
with respect to the M functions (pai^) 



12| and corrects the results of 
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and we find the they obey the following M coupled nonlinear differential equations (here we 
used the Lagrange multiplier e'^ to enforce Eq. (12. 2p ) 

E Ky€iy) = (2.9) 

with 

Ky = +\{5y,,+i-5y,,-,)+m{-lY5,y-\v:y, (2.10) 



xy 



1 



2ME^afc(2/'^) E 0r(^)0r(?/)j, (2.11) 

2iTi{x-y) ( a-b , ,\ 
1 ' e Ls K M ^"V 

Since Kah explicitly depends on 0, solving Eq. (12.91) is a self-consistent process which we 
perform by beginning with a random choice for p°^y, inserting it into h°j.y, diagonalizing the 
latter M matrices to obtain the M sets of eigenfunctions 0a(^)? ^'^d these then provide us 
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with a new estimate of p^^. This process is iterated until convergence is observed. The only 
restriction our initial random choice of p needs to obey is 

Y.Pl. = B + Lj2, (2.13) 

X 

which reflects Eq. (12. 3p and Eq. (12. 7p . Within the space of all functions that obey Eq. (12. 9p 
the correct solution is the one that has the lowest value of Tip. The latter is not equal to 
Yjan^ni siuce this will count the Coulomb interaction twice. Instead one has 

B+L s/2 M / \ 

2M 



Hsoiution/iv = ^ E E ( + E (c(^)0r + 1)) + m (-1)- 0:(x)0r (^)] 

n=l a=l \ X / 



4 i..M^ir4sin^(g(^ + 0/2)- ' 

Finally, to calculate the values of quark bilinears we build a single Dirac fermion from the 
staggered fields and what we get is 

M 

2M 



= -^E (pS,. + p"-m-i), (2.15) 



a=l 
M 



mix) = E (p:,x. - P^i,.-i) , (2.16) 



2M , 

a=l 
M 

M 
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(^^75^)(X) = E Imag {pl,^,) , (2.17) 

a=l 
1 M 

(V^7i V^) (X) = Real {pl,^,) , (2.18) 

a=l 

where here X = ax/2, x = even. In some of the plots that appear in forthcoming sections we 
present UV- regularized results for {ip'^'-)^ip){X). The regularization we choose is to subtract 
the corresponding values in the free theory. Specifically, anticipating large- volume inde- 
pendence, we subtract the free theory expectation values at infinite volume. In Appendix \K\ 
we show that this means we should perform the following subtractions 

(z^V)(X) ^ (^V)(X)-1, (2.19) 

(V'(1,Z75)^)(X) ^ (V5(l,Z75)^)(X)-(^(l,Z75)V')free, (2-20) 

- dv 1 

(^^)free = -m f- , (2.21) 

JpF ^TT ^ sin''(p/2) 

W'n.«.„ = +r'i , ""''"f ■ (2.22) 

Jpp ATT Jrn^ sin^(p/2) 



with pf = TT-j-^. Also, we always find that (tpjitp) = which means that time-reversal 
symmetry is intact. This, of course, also happens in the free theory (see the appendix).^ 

In the following sections we implement the numerical solution of the self-consistent 
Eqs. f l2A0l) -f l2:T2l) for different values of B, m/VA, L^A, and a^A. 

III. THE CASE OF B = AND B = 1. 

In this section we wish to see how sensitive are the values of the quark condensate and 
baryon mass to the technical parameter M. In jl^ we showed that the erroneous M = 1 
choice corresponds to setting the spatial Polyakov loops to unity while M = oo sets them 
to zero. Thus, we expect that the volume dependence in the former case to be strong, while 
in the latter it should disappear (for the connection between the expectation values of the 

u n 

Polyakov loops and volume dependence see [9|, [lOj. This expectation is confirmed by the 
results presented in Fig. [T] both the plots on the upper and lower panels of the figure were 



obtained for B = and a relatively fine lattice spacing of ayA/(27r) = 0.123. The plots 
present the way the UV-regularized value of {tptp)/N depends on the quark mass for different 
values of L. The only parameter by which these plots differ is M: it is fixed to 1 on the 
top panel and to 25 on the bottom. Fortunately, for the B = system there is an analytic 
prediction for arbitrary mass which we represent by the solid (red) curve in the figures 
and compare to our plots. As anticipated, while the L-dependence of the upper plot is very 
strong, it is hardly visible on the lower panel. This is large- Eguchi-Kawai independence 
at work. 

A glance at the lower panel of Fig. [T] might be alarming: our data seems to be quite 
far from the analytic solution in the regime of large quark masses. This, however, is a 
lattice artifact as we demonstrate in Fig. [2l There, we plot the m-dependence of the quark 



condensate for different lattice spacings, and for LJX/(2tt) = 32 (M was set to unity here 



since the 1/M corrections should be small for this volume). Indeed, we find that as ay A/(27r) 
drops from 0.123 to 0.0615 to 0.03075, the deviations of the data from the prediction of 



|15l | decrease. We take the deviations of the ayA/(27r) = 0.03075 data from the analytic 
prediction to reflect the numerical convergence error of our calculation, the 1/M, and 1/L 



The way we choose to build the continuum quantities is shghtly different from that of Ref. The 
difference, however, vanishes in the continuum hmit. 
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corrections, as well as the numerical error involved in evaluating the analytic expressions of 



Ref. 



15] . To check that these deviations are under control we plot in Fig. [3] the continuum 



extrapolation of (tptp) /{N^/X) for very light quarks with m/^/X = 0.05/v^27r, a relatively 



large volume of L^\/{2tx) = 32, and for several values of M. As the plot shows, while the 
data obtained for M = 15, 25 extrapolates to the analytic result in the continuum, the data 
of M = 1 does not. This discrepancy of ~ 0(8%) is thus a measure of the 1/M and 1/L 
corrections. 

Next we wish to examine the single baryon ground state. We begin by plotting, in 
Fig. m the baryon density {ip'^ ip) / N as a function of the spatial coordinate We choose 



a moderately large volume of L^X/ (27r) = 16, a lattice spacing of X/ (27r) = 0.123, and 



in Ref. 



present results for m/ J X/{2tt) = 0.05 and 1. These parameters are the same as those used 



ll|. Our calculation differs from the one in that paper by our choice of M : we 



use M = 25 while Ref. [ll| used M = 1. As seen in the Figure, the baryon becomes less 
localized as the quark mass drops. Eventually, at m = 0, it will spread throughout space, 
and become a massless delocalized objects. We see that the other quark bilinears, (ipip) and 
{ipi'j^ip): also modulate in space, but we postpone discussing this modulation to the next 
section. 

Also plotted in Fig. H] is the prediction for the soliton density of a single soliton in the 
sine-Gordon model. As argued in 3, [ll| this prediction should agree with the 't Hooft 
model's baryon density, but only for sufficiently small values of m/i/A. To emphasize the 



importance o: 
made in Ref. 



' having m/-\/A ^ 1 for this argument we mention that, in the approximation 
7, [ll|, the width of this soliton is determined by the value of (ipip) in the chiral 
limit, and one neglects all finite-m contributions to the latter. Remarkably, however, we see 
that the agreement between the prediction and our data is good even for m/ ^/X ~ 0.4, where 
the value of (ipip) differs by a factor of 4 from its value in the chiral limit. This agreement is 
not expected. For even smaller quark masses of m/^yX ~ 0.02 — 0.1 it was also a surprise to 



the authors of Ref. ll| . Here, and throughout the current paper, we see that such agreement 
seems good even for much heavier quarks than these. To make this point stronger we plot 
in Fig. Othe single baryon density for very heavy quarks of m/^/X ~ 1.2, where deviations 



Here, in contrast to Eq. (|2.4p , we denote by x the dimensional spatial coordinate (and not the lattice site 
index), i.e. from here on we denote the X of Eqs. (j2.15p - (|2.18p by x. 
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from the soliton density are already visible, but still very modest. Further surprising results, 
for even heavier quarks are presented in Section IIVI where we discuss finite density. 

To analyze the lattice corrections in our calculation we plot, in Fig. [6l the baryon density 



oi B = 1, m/y/X = 0.05/v/27r, M = 15, L^A/(27r) = 32, and for three different lattice 
spacings. The black continuous curve is again the number density of the one-kink ground 
state of the sine-Gordon model. As we see from the figure this prediction works very well as 



long as one compares to results of sufficiently fine lattice spacing. The data for a w A/ {2n 



0.03075 can be hardly distinguished from that of ay A/ (27r) = 0.067. We checked the lattice 
corrections for larger values of m, and found them to be smaller. 

As is clear from Figs. HHEl the single baryon ground state completely breaks translation 
invariance. That this symmetry can break is a result of taking the large- limit : the energy 
cost of moving the baryon is proportional to A^, and the latter is strictly infinite in our work. 
For finite values of A^ the true ground state will presumably have total zero momentum and 
will not break translations. A finite density of baryons can, however, induce a spontaneous 
breakdown of translations, even at finite A^. This will happen if the ground state contains a 
crystal of baryons which leaves invariant only a subgroup of the full translation symmetry. 
The cost of moving such a crystal scales like (A^ x Volume) and so will be infinite in the 
thermodynamical limit, for finite values of A^ as well.^ ^ In any event, we work at A^ = cxd, 
and so in our case even a single baryon breaks translation invariance. 

We proceed by analyzing the baryon mass rriB which we calculate by subtracting the 
energy of the B = grou nd state from the energy of the B = 1 ground state. Again, this 



was first done in Ref. HI], but by using M = 1. Here we wish to check whether the volume 



chosen there was large enough for the 1/M corrections to be negligible. The results are 



presented in Fig. [7] where we see that for LyA/(27r) = 16 and ayA/(27r) = 0.123, the 1/M 
corrections to itlb/ {N\/~X) are of ~ 0(5%) for a light quark mass of m/^fX ~ 0.02, and 
around ~ 0(1 % — 2 %) for heavy quark masses of m/ \/~\ ~ 0.4. As a comparison, for the 
parameters chosen in the lower panel of Fig. [3, Ref. [11] obtained niB/iNy/X) ~ 0.579, while 
we see that the large-M extrapolation results in mB/{N\/X) ~ 0.592. 



^ We thank Yigal Shamir for emphasizing this issue to us. 

^ In our one-dimcnsional case, any spontaneous breaking of continuous symmetry is removed by IR fluctu- 
ations, but these are suppressed as long as we are at large- A^. 
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A. Further comparisons with the sine-Gordon soHton 



As mentioned above, for light quarks obeying m/^fX <^ 1 the continuum version of 
the classical Hamiltonian Tip reduces, for a certain ansatz for ,, to the Hamiltonian of 

|. For B = I this ansatz is the sine-Gordon soliton and in this 
subsection we wish to check our calculations further by comparing the energy and height of 
the sine-Gordon soliton with itlb and with the maximum of the baryon density Pmax- 

We begin by performing continuum extrapolations for and Pmax obtained for three 
values of the quark mass : m /y/X = (0.5, 0.25, 0.05)/v^. These results were calculated 



using M = 15 and L^\/{2'k) = 16,24,48 respectively, and are presented in Figs. [H]and[9l 
In Fig. [To] we compare the continuum values of and Pmax to the following predictions of 
the sine-Gordon model. 



( Y _ 8 m / (^V)max y _ 2 m 

As the figure shows this prediction agrees with our data for small enough masses. Specifically, 
this means that the baryon mass and the inverse of its width, decrease with the quark mass 
to zero, and that that they so approximately like the Sine-Gordon model predicts, i.e. linear 
in y/m. Thus, at small masses, the baryon turns into a delocalized massless excitation of 
the vacuum. This is very different than the way the four-dimensional baryon behaves in the 
chiral limit, and is precisely the reason we wanted to solve the 't Hooft model for general 
quark masses. 

IV. FINITE DENSITY AND THE EQUATION OF STATE 

This section contains the main results of our study: the ground state of a system at 
nonzero baryon density. The densities we study are obtained by putting a large baryon 
number of < 5 < 30 in a box of fixed size L. We begin by presenting, in Figs. [TTHTHl 
the results of the quark condensate {il)il))/N and the density {il)'^il))/N as a function of the 
spatial coordinate for different values of the quark mass, the lattice spacing, the volume, and 
B. What we find is that for any nonzero value of 5, translation symmetry is spontaneously 
broken and that at large enough values of B/L, a crystal is formed. 

The crystal that we find is a direct result of the repulsion between the baryons. This 



repulsion was mentioned as an aside in Ref. 



ll|, and in this paper we show it explicitly when 



we discuss the equation of state (see the next subsection). In higher dimensions one also 
expects that a formation of a crystal, as long as is infinite To understand this think 
about the baryons as classical interacting particles (rather than emergent coherent objects 
that do not fluctuate). In this picture a baryon-baryon potential can either be repulsive 
or attractive, depending on the dynamics. In the former case (which is what happens in 
the single flavor 1 + 1 case, as we show in our paper), putting a bunch of baryons in a box 
obviously creates a crystal. What happens in four dimensions? physical, 3-color, baryons 
have long-distance attraction in four dimensions, but there is still a hard-core repulsion at 
short distances. Thus, the distance at which the baryon-baryon potential has a minimum, 
is a finite, O(Aqcd) value, which we denote by tq. If this picture survives the large- limit, 
then the absence of fluctuations at that limit tells us that the inter-baryon distance will 
exactly be tq, which means that they will form a crystal at infinite volume. 

Naturally, the crystal structure is contaminated with lattice artifacts when its wave length 
is close to the lattice spacing (see the odd looking structure of the crystal for i? = 30 in the 
right panel of Fig. [TTl) . These artifacts go away when we either decrease the density (see 
S = 10 in the same figure) or decrease the lattice spacing (see 5 = 30 in the right panel of 
Fig. [T2|) . This phenomenon, of increasing lattice artifacts with increasing density, is of course 
anticipated in advance, and is usually referred to as 'lattice saturation': it happens when the 
number of baryons per site is close to 1, which is the maximal number that Pauli-exclusion 
principle allows for a theory with a single flavor. 

In Fig. [14] we plot {ip i'^^ip) /N, and find that it is also spatially modulated. In Fig. [15] 
we show that this crystal has a helical structure: both {%l)il))/N and {tp ij^ip) / N are plotted 
versus the spatial coordinate. The parameters of the data plotted are the same as those used 
to generate Fig[Tl] In Fig. [TH] we repeat the plot presented in Fig. [TJl but for a lighter quark 
mass of m/^/X = 0.05/-\/27r. In Fig. [T7]we present the baryon density for the light quark 
case and compare to the predictions of the sine-Gordon crystal - see [7|. The agreement is 
impressive. We present similar plots for heavier quark masses of m/v^~ 0.2,0.4, 1.2,2.4 m 
Fig. [IHJ The deviations from the sine-Gordon crystal are increasing with increasing m/i/A, 
but at a relatively modest rate. 

Next, we integrate {ijjip) over the volume and present the results versus the density in 
Fig. [Tn]-[2ni What we find is that the volume averaged value of (ipip) drops with the density 
in a way that becomes stronger with decreasing quark mass. 
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A. The equation of state 

Moving on, we plot the energy E of the ground state in the presence of B baryons, i.e. 
the equation of state, in the upper panel of Fig. [2T1 Here we regularized E by subtracting 
the B = ground state energy. We present results for both light and moderately heavy 
quarks (see captions). Also plotted in Fig. [211 as a solid (black) curve, is the equation of 
state for free massless quarks, in the continuum. The latter is given by 



and seems to be quite close to our data, especially for large enough density. This means 
that there is a smooth transition from Baryon physics to free quark physics as the density 
increases.^ This seems to be a nontrivial result and it would be very interesting to understand 
it better. 

Note that the fact that the E{B) is quadratic in B at large- -B means that the would-be 
binding energy A, which is given by 



is proportional to B, is growing, and is always positive. This means that there is no binding 
in 1 + 1 and a single flavor, or put differently that the baryons repel each other. This also 
means that they would create a crystal when put together in a box, as we showed in the 
previous subsection. 

In the lower panel of Fig. [2T]we present the numeric evaluation of dE{B)/dB which is 
equal to the chemical potential. Here, in contrast to the upper panel, the four data sets 
do not fall on top of another because of the different volumes used. Note that we have 
normalized the y-axis in this plot to the constituent quark mass itib/N . 

Finally, in Fig. [22] we combine Figs. [IHl [20], and [21], and present the values of the regu- 
larized quark condensate, divided by its 5 = value, and plotted vs the chemical poten- 
tial (the latter divided by m^/A^). The absence of data points between fi/{mB/N) = 
and ^/ {itlb/N) = 1 reflects the absence of excitations with nonzero baryon number that 
are lighter than a single baryon and thus the /x-independence of the ground state for 
0< I2< tub/N. 

^ We thank D. Son for pointing this issue to us. 




(4.1) 




(4.2) 
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V. PARTIAL VOLUME INDEPENDENCE AT LARGE- 



T he g oal of this section is to check the following simple observation that was made in 
Ref. 1(|: Despite the breakdown of translation invariance in the presence of baryons, there 
is a remnant of the "Eguchi-Kawai" volume independence of large-A^ QCD which allows 
us to calculate properties of B baryons in a box of size L by putting a single baryon in a 
box of size L/B. To see this 'soft' form of large- volume independence the spatial gluonic 
zero modes need to be treated correctly, and as explained above this amounts to using large 
values of M. To check this observation we compare our numerical solution of a system with 



B = 5 baryons and a volume of LyA/(27r) = 16 to the B = 1 case and L' = L/5. We 
present the results for the quark condensate and density obtained with M = 1 in Fig. [23] 
and with M = 25 in Fig. [211 As the plots clearly demonstrate such volume independence 
takes place only for large values of M (i.e. as long as one treats the zero modes correctly). 

Clearly, the observation we make in this section can have an important practical con- 
sequence: if it works in four dimensions as well, it can save computational resources for 
lattice Monte-Carlo simulations. To see this observe that what determines the size of the 
lattice corrections is the ratio a/ A, where a is the lattice spacing and A is the crystal wave 
length. Thus, to put more baryons in a box, we need to decrease the lattice spacing, and 
fixing the physical volume, this means increasing the number of lattice points, and thus the 
calculational cost. If however, we assume that the physical baryon is close to the infinite- iV 
baryon, then according to what we find in this section, we can study large densities by 
putting a single baryon in a box whose physical size is small. This can be accomplished 
by just decreasing the lattice spacing, but by keeping the number of lattice sites fixed. At 
any given volume, the ratio a/ A is equal to a/L, but the latter is equal to l/Ng, with A^^ 
denoting the number of sites in one of the lattice directions, which is kept fixed and large. 

The usefulness of this proposal is predicated on two assumptions. First we assume that 
physical, 3-color, QCD is well approximated by its large-A^ limit, even in the presence of 
baryons. We do not see any obvious problems with this assumption. Second, it assumes that 
in four dimensions the ground state of the dense system is a crystal. Moreover, it assumes 
that the crystal structure is commensurate with the lattice topology (i.e. that it is a simple 
cubic). As we say in the previous section, because we are at large- a crystal is expected 
to form in any number of dimensions, and at finite-A^ this crystal may dissolve, but can 
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leave behind a nontrivial structure in correlation functions, that can be studied with the 
suggestion we make in this section. 

The question of whether the formed crystal is simple-cubic is a dynamical one and remains 
an unproven assumption in our context. 

VI. SUMMARY AND CONCLUSIONS 

In this paper we studied large-A^ QCD in 1 + 1 dimensions and solved for its ground 
state given arbitrary baryon number B, quark mass m, and spatial volume L. We used 
the Hamiltonian lattice formalism and regularized the theory in the IR by placing it in a 
finite spatial box. This IR regularization prevents us from gauging away the spatial gluons 
completely, and a set of gluonic zero modes remains and plays a crucial role at small volumes 
and/or large baryon densities. These zero modes were ignored in some studies of this theory 
at i? = and B = 1 and here we analyzed the effect such a mistreatment has on the value 
of quark condensate (ipip) at B = 0, and on the baryon mass. Next we analyzed the ground 
state of a variety of baryon numbers in the range < 5 < 30, and how it depends on the 
volume, on the lattice spacing, and on the quark mass. We find that these ground states 
break translation invariance and form a crystal for large enough densities, and that this 
crystal is characterized by a helical structure in the quark bilinear condensate densities. For 
light quarks our results closely resemble those obtained analytically by Schon and Thies 
in Ref. [3| who used an ansatz that reduces the 't Hooft equations, for mj ^\ <^ 1, to 
the corresponding equations in the sine-Gordon model. This work also ignores the gluonic 
zero modes (which seems justifiable since it focused on the infinite volume limit). In any 
event, such an agreement is not expected a priori because, first, we make no ansatz for the 
ground state (but rather minimize the classical Hamiltonian of the system numerically), 
and, second, most of the quark masses we explore are not light, and some are quite heavy 
with m/^yX ~ 1.2 - 2.4. 

In general, and for any value of m > 0, our results show that the ground state of large- 
QCD in 1 + 1 dimensions strongly depends on the baryon number B. In particular, for 
any B > 1 the ground state breaks translations. To see how this dependence is refiected 
in the grand-canonical ensemble, we used the way the ground state energies depend on the 
baryon number and calculated the quark chemical potential /i. This results in the phase 
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structure sketched in Fig. [22] : there is a phase transition that occurs at a value of n equal to 
the baryon mass (divided by A^) and that separates a phase which has zero density, that is 
//-independent, and that is invariant to translation symmetry, from a nonzero density phase 
with spontaneously broken translation symmetry and in which physical observables strongly 
depend on /i. 

Finally, it is interesting that when we calculate the equation of state, we find that for 
large-i?, the energy grows quadratically in B. This also means that the baryons repel each 
other, which explains the formation of the crystal. We also see that this quadratic growth 
is quantitatively quite close to the quadratic growth in the equation of state of free massless 
quarks (see Fig. [2T] and Section TlV Al) . Thus it seems that for a single flavor and in 1 + 1 
dimensions there is a smooth connection between baryon physics at small/moderate density 
and free quark physics at large densities. 



As mentioned in the introduction, Ref. 



has attempted to approach large-A^ two- 



dimensional QCD by studying, instead of 't Hooft model itself, its associated 'zero-volume 



Eguchi-Kawai (EK) matrix model {9!]. The main conclusion of Ref. js] is very puzzling: 
large-A^ QCD in 1 + 1 dimensions is argued to be completely independent of the baryon 
chemical potential. This, if correct, contradicts simple physical intuition (and of course 
contradicts all the results of Refs. [ill and of the current paper). Specifically, what the 
conclusions of Ref. mean is that there are no baryons in two dimensions and so this also 

n 

contradicts section IX of Ref. pL7| • 

One of the reasons why the conclusion of Isj is questionable is the following. The use of 
the EK prescription by the authors of Ref. [8| assumes that the euclidean Dirac operator 
factorizes in momentum space. This assumption is predicated on having a ground state 
that is invariant under translations. If the latter symmetry is spontaneously broken in 
the ground state, then a quark can change its momentum by interacting with momentum- 
carrying condensates, and such interactions invalidate the factorization property assumed by 
the authors of Ref. [8|.^ Put differently, one can easily show that the generating functional of 



momentum-carrying quark bilinears, such as / dx {ip{x) ip{x)) e**^^, which 



cannot be calculated with the momentum factorization used by Ref. jsj. Indeed, general 



Dreak translations. 



That the Dirac operator does not factorize in momentum space when translations are spontaneously 
broken is well-known and we refer the reader to [3] for explicit examples. 
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arguments tell us that Eguclii-Kawai reduction will take place only so long as the ground 
state is translation invariant (see for example the first paper of Ref. [3]), and since we show 
that for fi > rriB/N translation symmetry spontaneously breaks, then one cannot use the 
EK matrix model to study the field theory in that regime. 

n □ 

Another reason why one may question the conclusions of Ref. [8[ was exposed in Ref. [10| 
and is related to general features of the way baryon density would behave within single-site 
models. The issue is simply stated: if we have any nonzero baryon number i? > 1 on a 
lattice that has a single site, then, by construction, we force the baryon density to be at the 
cutoff scale, i.e. of 0(l/(lattice — spacing)). This is true for any lattice coupling, and does 
not get 'better' as we take the naive continuum limit by sending the bare lattice spacing to 
zero. Thus, by their definition, single-site models can either have zero density, or density 
at the cutoff scale, and no densities in between. Indeed, this explanation fits well with the 
findings of Ref. 8j, where a transition from a zero density phase to a nonzero density phase 
with 0(l/(lattice — spacing)) density is reported. This was also seen in [19']. While Ref. js] 
interprets this as having no dependence on the chemical potential in two-dimensional QCD, 
we suggest that by working on a single site, the approach of Ref. js] is 'forcing' the system 
to show only this unphysical behavior. It seems reasonable that the real, physical, transition 
would be seen only if one had a finite number of sites Lg, and that to keep lattice artifacts 
under control, one would have to send Lg to infinity. With such a finite number of sites, one 
may also be able to see the breakdown of translation symmetry. 
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APPENDIX A: FREE THEORY SUBTRACTIONS 



For the purpose of subtracting the UV divergences from expectation values, we calculate 
the values of the quark condensate for the free theory in this appendix. Our starting point 
is the differential equation for in the free theory, which we obtain by setting A = into 
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Eq. (ETTUj) . This gives 



- + 1) - - 1)) + m (-ir 



(Al) 



It is convenient to define a new lattice with double the lattice spacing and a unit cell that 
contains two fields. We choose the following convention (from here on we suppress the index 
a since it enters Eq. (lAll) trivially) 



X = even, 
X = odd. 



(A2) 



In terms of (f)e,o{X), Eq. fl2.10l) becomes the following two sets of coupled equations 

mx) - - 1)) - mc(x) = 

+ 1) - 0^(X)) + m0:(X) = eirAX). 

Since X obtains the values 1, 2, 3, ... , Ls/2 we Fourier transform it as 

1 



(A3) 
(A4) 



(A5) 



^iv,/2 p 

where here p = 27rfc/(L^/2) and A; = 1, 2, 3, ... , Ls/2. Substituting this into Eqs. flA3ll-flA4l) 
gives the following matrix equation 



-m 



_g+ip/2 sin(p/2) 




^0^(p) 



(A6) 



whose orthonormal eigenvectors are 



_p-ip/2 



sin(p/2) 



2v/m2 + sin2(p/2) ( Jm^ + sin2(p/2) ± m] V ^ ± V^' + sin (p/2) 



(A7) 



and whose eigenvalues are = iym^ + sin^(p/2). 

For a fixed value for the baryon number the elements of the matrix density pxy are set 
by Eq. (12.71) and in our case the index n in that equation corresponds to the combination of 
the momentum index p and the ± index of Eq. (1A7I) . Since we are seeking for the solution 
with minimum energy, the sum over the combined index n = (p, ±) that forms p^y (see 
Eq. (12.71) ) includes the following : 



16 



• All the eigenvectors (p) with p = 1,2, . . . , Ls/2 (all these have negative energy e (p) 
and so by including all of them we lower the energy to a minimum) . 

• A set of B eigenvectors (f)'^{p) whose energy e~^{p) is the lowest. Below we denote the 
set of momenta that corresponds to these eigenvectors by F. 

Focusing only on the diagonal elements of p, this results in the following form 



Ls/2 



ElC(p)P + E \<Pt{p)\' 

P p&F 



X = even, 



X = odd. 



(A8) 



Using Eqs. fl2.15p -( l2.16l) . and after some algebra, we obtain 

(^V)free(X) 



r E+E =5/^^ + ^, 



(^^)f,ee(X) 



m 

1 

Ts 
1 



E 



1 



E 



1 



E 



p \im? + siv?{p/2) pgF ym^ + sin^(p/2) 
sin^(p/2) ^ sin^(p/2) 

rh? + sin^(p/2) peF Jm? + sin^(p/2) 



(^7l^)free(X) = ^ ^ E 



sin(p) 



+ E 



sin(p) 



(A9) 
(AlO) 
(All) 
(A12) 



p \lm? + sin^(p/2) pgF ^Jrh? + sin^(j9/2) 

Eq. ( 1A9P tells us to subtract 1/2 from the density of the interacting theory, hence Eq. (12.191) . 
In the infinite volume limit the remainder is 



(VV(^))free = \ f_" ^ = B/L, 



which gives us the value of the Fermi momentum pp = tt 



B 
Ls/2- 



(A13) 



Using this and taking the 



infinite volume limit of Eqs. [XToMTTl we obtain Eq. (IXZTD and Eq. (12:221) . Finally, the 
symmetry p —p of the sums in Eq. (IA12p tells us that {ijj 71 ip) = 0. 
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FIG. 1: Quark mass dependence of the condensate for B = and different volumes. Upper panel: 
Mistreating the spatial gluonic zero modes (setting M =1). Lower panel: Incorporating the zero 
modes with M = 25. The red curve is the result of Ref. [15| (see text). As explained in the text, 
the discrepancy between our data in the left panel and the analytic curve reflects lattice artifacts. 
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FIG. 2: Quark mass dependence of the condensate for B = Q and different lattice spacings. We 



use M = 1, but a relatively large volume of L^J\/{2ti) = 32, where the zero modes do not 
induce any significant 1/M effect. The red curve is the exact result in the continuum and infinite 



volume limit from 



151], which agrees very well with the results from the smallest lattice spacing of 



a^\/{2TT) = 0.03075. 
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FIG. 3: Continuum extrapolation of the quark condensate for mlyX = 0.05/\/27r and for several 

n 

values of the parameter M. The Filled circle (black) at a = is the analytic result [15[. All data 
was obtained for the same volume of L\/ A/ (27r) = 32. 
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Baryon density of a single baryon 



□ mNX=M^{2%) 

o mNX=0.05N{2K) 

sine-Gordon soliton for mNX = 1/V(27i) 

sine-Gordon soliton for mNX = 0.05/V(27r) 
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FIG. 4: Baryon density of a single baryon. Squares (blue): an intermediate quark mass; m/^/X = 
l/\/27r. Circle (green): Small quark mass; mj^fX = 0.05/\/27r. The dashed (blue) and solid (green) 
curves are the predictions of the sine-Gordon soliton that are expected to be descriptive of the data 
in the m/^A limit. 
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Baryon density of the B=1 ground state, mHX ~ 1 .2 



0.45 




FIG. 5: Baryon density of a single baryon obtained for very heavy quarks : m/^/X = 3/-v/27r. 
Again the curve is the prediction of the sine-Gordon sohton vahd in the hght quarks regime of 
m/^A <C 1. 
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Baryon density of a single baryon 
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FIG. 6: Baryon density of a single baryon for different lattice spacings. The solid curve (black) is 
the result of the sine-Gordon. Note the different scales of the x and y axes compared to those of 
Fig.H 
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FIG. 7: Baryon mass at a^/X/(2^T) = 0.123 and L^J \/ [2ii) = 16 as a function of 1/M. 
Top panel: An intermediate quark mass; m/\/A = l/\/27r. Bottom panel: Small quark mass; 
m/\/A = 0.05/^2^. 
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FIG. 9: Same as Fig.[8]but for (Max(V't^)) /{N^X). 
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FIG. 10: Comparison of the results obtained from the continuum extrapolation appearing in Fig. [8] 
and Fig. [9] to the predictions of the sine-Gordon soliton. 



Quark condensate, m HX ~ 0.2 Number density, m NX ~ 0.2 




xVA, x^X 



FIG. 11: Expectation values of {■iIjiI^)/N (left panel) and the density {il)^ip)/N (right panel) as a 
function of the spatial coordinate x. Here ay^ A/(27r) = 0.123, m/^/X = 0.5/-v/27r, A/ {2tt) = 16, 
and M = 25. The baryon number showed are B = 0,3, 5, 10, 15, 20, 25, 30. 
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FIG. 13: Same as Fig. ([II]), but for a^X/{27r) = 0.246, m/^/X = 0.05/^/2^, L^/X/{2nj = 32, and 
M = 15. Note the different scale in the x-axis compared to the one of Figs. [TTHT2] 
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FIG. 14: The condensate (ipij^ip) /Nfor B = 4 and an intermediate quark mass of m/^/X = 
0.5/V27r. To produce this plot we used aVA/(27r) = 0.123, L^A/(27r) = 16, and M = 25. 
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Helical structure of the B=4 ground state, m /V^ ~ 0.2 

a ' q 
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(Chiral condensate) / x-^X 



FIG. 15: The helical structure of the ground state of S = 4 and an intermediate quark mass of 
m/\/A = 0.5/\/27r; (ipil^) on the y-axis, (ipij^il^) on the z axis, and the spatial coordinate x\/A on 
the X axis. Technical parameters are the same as those of Fig. [TH 
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Helical structure of the B=4 ground state, m^NX ~ 0.02 
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(Chiral condensate) / xV?i 



FIG. 16: Same as Fig. (jlSp . but for rn-/\/A = 0.05/\/27r. Other technical parameters are 
a^/X/(2^F) = 0.123, L^J\/{2ti) = 32, and M = 15. Note the different scales of the x, y and z 
axes compared to those of Fig. [151 
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Baryon density of the B=4 ground state, m^NX ~ 0.02 
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FIG. 17: The baryon density for m/^/X = 0m/^/2^, B = A, Ly/X/{2n) = 32 and a^A/(27r) = 
0.123. The sohd (black) line is the prediction of the sine-Gordon crystal. The technical parameter 
M was fixed to 15. 
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FIG. 18: The baryon density for B = A. Upper-left panel: intermediately heavy quarks with 
m/\/A = 0.5/\/27r and A/ (27r) = 0.0615. Upper-right panel: heavy quarks with m/\/A = 
1.0/\/27r and a\/A/(27r) = 0.123. Lower-left panel: heavier quarks with m/\/A = 3.0/\/27r 
and a-\/ A/(27r) = 0.075. Lower-right panel: even heavier quarks with m/\/A = 6.0/\/27r and 



a-\/A/(27r) = 0.075. The solid (black) lines are the prediction of the sine-Gordon crystal. The 
parameter M was fixed to 25 for the upper-left panel, and 15 for the other panels. In units of 



^A/(27r), the volumes were 16 for the upper plots, 10 for the lower-left plot and 8 for the lower-right 
plot. 
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Volume averaged Quark condensate 




B / (LV^i) 

FIG. 19: The volume average of ('0V') vs the baryon density B/L, for a^J\/{2Ti) = 0.123 and 
L^/XJ(2tt) = 16; Circles are for ?Ti/\/A = 0.5/\/27r and squares are for m/^/X = 0.05/V2vr. For 
each mass we show results from both M = 1 and M = 15 (open red and closed blue symbols). 
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Volume averaged Quark condensate, m NX ~ 0.2 
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FIG. 20: The volume average of {ipip) for M = 1 (open symbols) and M = 15 (filled symbols) and 
L\/\/{2t^ = 16; Circles are for a^J A/ (27r) = 0.123 and squares are for a^J A/ (27r) = 0.0615. Here 
the mass is m /VA = 0.5/^/2^. 
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FIG. 21: Upper panel: The S-dependence of the difference in energy AE{B) between the B > and 
the B = systems. Lower panel: The chemical potential normalized to the baryon mass divided 
by N. For these calculations we choose M = 15 and the other details of the calculations appear 
in the legend. The value of the baryon mass for the light (m/\/A = 0.05/\/27r)and intermediate 
(m/\/A = 0.5/\/27r) quark masses is niB/ {'sfXN) ~ 0.11,0.37, respectively. 
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Volume average of quark condensate vs chemical potential 
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FIG. 22: The value of the quark condensate for m/^/X = 0.5/\/27r, L^A = \Q^J2'K, and a^J\ = 
0.123V27r (green squares) and m/y/\ = 0.05/^/2^, L^A = 32^27r, and a^X = 0.246V27r (blue 
circles). The x-axis is the chemical potential (normalized to the baryon mass), and the y-axis is 
the regularized condensate, normalized to its value in the B = Q system. 
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0.15 



Quark condensate, m NX ~ 0.2 
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Number density, m^/VX ~ 0.2 




FIG. 23: The way (^/'V') (top panel) and the baryon density (bottom panel) depend on the spatial 
coordinate for M = 1 (i.e. by mistreating the zero modes). Results in dots (blue) are for a system 
of S = 5 in a volume of L^/\/{2tt) = 16 and results in squares (red) are for B = \ and L' = L/5. 
Clearly the mistreatment of the zero modes, reflected by the choice M = 1, does not lead to volume 
independence - the squares and the dots do not coincide. 



38 



0.15 



Quark condensate, m NX ~ 0.2 
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FIG. 24: The same as Fig. (j'23|) . but when one correctly incorporates the zero modes in the 
calculation : here M = 25 rather than M = 1. In this case the dots and squares coincide and 
partial large-A^ volume independence holds. 
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FIG. 25: The phase diagram of two-dimensional large-A^ QCD in the grand-canonical ensemble. 
For quark chemical potentials fi that obey fi > ms/N, translation symmetry is broken and the 
ground state depends on fi. 
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